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INTRODUCTION 


When  numerically  simulating  an  event  such  as  the  penetration  of  an  armor,  where  mate¬ 
rial  deforms  at  a  high  strain  rate  and  temperature,  results  will  significantly  depend  upon  the 
choice  and  the  validity  of  the  material  constitutive  model  that  is  used  in  the  simulation.  To 
develop  new  models,  evaluate  existing  models,  and  to  obtain  model  parameter  constants,  it  is 
necessary  to  obtain  experimental  stress-strain  data  as  high  deformation  rates. 

One  method  of  obtaining  experimental  high  strain  rate  (100  to  104  s'1 2 3 4 5 6 7 8 9),  stress-strain  data 
is  by  using  the  split  Hopkinson  (also  known  as  Kolsky  bar)  bar  technique.  Split-Hopkinson 
bar  configuration  consists  of  two  elastic  bars  between  which  the  specimen  is  held.  In  this 
test  system,  a  stress  pulse  is  applied  to  one  bar  (input  bar)  in  the  direction  of  the  other  bar 
(output  bar)  by  an  external  mean;  for  example,  by  propelling  a  striker  bar  at  the  input  bar  or 
releasing  a  stored  torque  in  the  input  bar.  Part  of  the  stress  pulse  will  travel  to  the  output 
bar  (through  the  specimen)  while  the  remainder  will  reflect  back  into  the  input  bar.  By 
measuring  the  elastic  stress  waves  in  the  input  (reflected)  and  output  (transmitted)  bars, 

Kolsky  showed  that  it  is  possible  to  determine  the  stress-strain  response  of  the  sandwiched 
specimen.1  It  has  been  demonstrated  that  this  technique  can  be  used  to  obtain  stress-strain 
response  under  compressive,  ’  tensile,  '  or  torsional  ’  conditions. 

Though  the  high  strain  rate  behavior  can  be  obtained  by  either  compression  or  tension, 
there  are  many  advantages  to  obtaining  the  stress-strain  behavior  by  the  torsional  (shear 
mode)  method.  For  example,  when  testing  for  large  strains,  the  strains  that  can  be  achieved 
by  tension  is  limited  by  necking,  and  compression  is  limited  by  barreling  of  the  specimen.  In 
addition,  errors  arising  from  dispersion  of  longitudinal  stress  wave  on  its  travel  along  the 
elastic  bars  in  a  compression  or  tension  test  are  not  present  in  a  torsion  test;  there  is  no 
wave  dispersion  in  shear  wave  transmission.  Although  it  is  assumed  that  the  stress  conditions 
are  uniaxial  in  a  compression  or  tension  test,  at  high  rates,  due  to  radial  inertia  (radial  stress) 
that  is  opposing  the  Poisson’s  expansion  associated  with  uniaxial  stress,  stress  conditions  may 
not  be  fully  uniaxial.  In  a  compression  test,  at  high  strains  and/or  strain  rates,  this  inertia 
effect  may  get  accentuated  by  the  radial  friction  at  the  specimen  and  bar  interface,  even  with 
the  presence  of  lubrication.  However,  in  spite  of  these  disadvantages  and  limitations,  the  com¬ 
pression  test  has  the  advantage  of  requiring  a  simpler  and  cheaper  specimen  compared  to  the 
torsion  or  tension  test. 

Shear  also  is  the  main  form  of  deformation  that  is  present  in  high  rate  deformation 
events  such  as  penetration  of  an  armor.  Therefore,  to  accurately  predict  material  behavior  in 
these  types  of  situations,  material  constitutive  models  should  be  developed  from  shear  data  at 


1.  KOLSKY,  H.  Stress  Waves  in  Solids.  Oxford  University  Press,  London,  1953. 

2.  LINDHOLM,  U.  S.  Some  Experiments  with  the  Split  Hopkinson  Pressure  Bar.  J.  Mech.  Phys.  Solids,  v.  12,  1964,  p.  317. 

3.  FRANTZ,  C.  E.,  FOLLANSBEE,  P.  S.,  and  WRIGHT,  C.  E.  New  Experimental  Techniques  With  The  Split-1  lopkinson  Pressure  Bar.  in 
High  Energy  Rate  Forming  -  1984,  Berman,  I.,  and  Schroeder.  Ed.,  ASME,  New  York,  1984,  p.  229. 

4.  FOLLANSBEE,  P.  S.  The  Hopkinson  Bar,  High  Strain  Rate  Compression  Testing,  in  High  Strain  Rate  Testing,  Mechanical  Testing.  Metal 
Hand  Book,  American  Society  for  Metals,  v.  8,  p.  198 

5.  LINDHOLM,  U.  S.,  and  YEAKLEY,  L.  M.  High  Strain  Rate  Testing:  Tension  and  Compression.  Exp.  Mech.,  v.  8,  1908.  p.  1. 

6.  HARDING,  J.,  WOOD,  E.  D.,  and  CAMPBELL,  J.  D.  Tensile  Testing  of  Material  at  Impact  Rate  of  Strain.  J.  Mech.  Eng.  Sci.,  v.  2.  19o(), 

p.  88. 

7.  NICHOLAS  T.,  and  BLESS,  S.  J.  Split-Hopkinson  Bar  in  Tension,  High  Strain  Rate  Tension  Testing,  in  High  Strain  Rate  Testing 
Mechanical  Testing.  Metal  Hand  Book,  American  Society  for  Metals,  v.  8,  p.  212. 

8.  BAKER,  W.  W.,  and  YEW,  C  H.  Strain  Rate  Effects  in  the  Propagation  of  Torsional  Plastic  Waves.  J.  Appl.  Mech.,  v.  33,  1966,  p.  917. 

9.  HARTLEY,  K.  A  .  DUFFY,  J.,  and  HAWLEY,  R.  H.  The  Torsional  Kolsky  (Split-Hopkinson)  Bar,  High  Strain  Rate  Testing  Mechanical 
Testing.  Metal  Hand  book,  American  Society  for  Metals,  v.  8,  p.  218. 
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high  strain  rate  and  temperature.  In  addition,  test  results  converted  to  shear  from  tension  or 
compression  data  by  using  a  criterion  such  as  von  Mises  equivalent  stress-strain  is  not  valid  at 
large  strains  above  approximately  20%. *'10  Shear  tests  can  be  conducted  using  a  punch  or 
double-notch  type  shear  tests  using  a  compression  bar.11  In  this  case,  although  high  strain 
rates  in  the  order  of  104  s'1  could  be  achieved,  disadvantages  with  compression  bar  such  as 
geometric  wave  dispersion  associated  with  longitudinal  wave  transmission,  as  well  as  non-uni- 
form  deformation  in  the  specimen  gage  section  will  be  present.  However,  a  torsion  test  elimi¬ 
nates  most  of  the  above  disadvantages.  It  is  noted  that  rate  change  tests,  also  known  as 
jump  tests,  (change  of  state  keeping  state  variables  representing  the  structure  constant)  that 
are  required  to  develop  and  evaluate  modern  physically  (dislocation  dynamics)  based  high  rate 
and  temperature  state  variable  constitutive  models  can  be  easily  conducted  using  the  torsional 
method.  In  a  forsion  test,  shear  stress  pulse  could  be  generated  and  released  by  explosives. 
This  method  of  torsion  testing  could  generate,  in  addition  to  the  torsional  pulse,  a  bending 
and  axial  pulse,  as  well  as  ringing  in  the  torsional  pulse,  thus  requiring  elaborate  mechanical 
filtering  methods.  In  addition,  an  explosive  capability  has  to  be  set  up  in  the  laboratory.  In 
contrast,  a  torsional  pulse  that  is  generated  and  released  by  mechanical  means  (with  mechani¬ 
cally  releasing  a  clamp  that  would  hold  a  stored  torque  that  is  generated  by  twisting  the  input 
bar)  will  not  have  the  above  complications,  but  will  have  higher  rise-time  and  lower  maximum 
achievable  strain  rate. 

To  facilitate  the  development  and  evaluation  of  constitutive  models  that  would  better  pre¬ 
dict  shear  behavior  at  high  strain  rate  and  temperature,  and  to  make  measurements  such  as 
strain  and  temperature  in  shear  bands  while  they  are  being  generated,  a  Stored  Torque  Tor¬ 
sional  Hopkinson  Bar  (THBar)  was  developed  at  the  U.S.  Army  Materials  Technology  Labora¬ 
tory  (MTL).  This  report  presents  the  information  and  characteristics  of  the  Torsional 
Hopkinson  bar  that  was  developed  at  MTL.  Some  preliminary  torsional  stress-strain  data,  at 
high  strain  rate  and  room  temperature,  for  annealed  OFHC  copper  (99.99%  copper)  from  this 
test  instrument  is  included  in  this  report.  High/low  temperature  and  jump  test  capability  arc 
being  added  to  the  THBar. 


DESCRIPTION  OF  THE  APPARATUS 

The  design  principal  of  our  THBar  is  similar  to  that  of  the  machine  located  at  Brown 
University.9  Here,  the  THBar  consists  of  two  elastic  bars,  a  clamp  with  a  sudden  releasing 
mechanism,  and  a  loading  arrangement  to  twist  one  of  the  bars.  Each  end  of  the  specimen 
is  attached  to  one  end  of  each  bar  so  that  the  specimen  is  held  between  them.  The  clamp  is 
used  to  hold  the  bar  from  rotating,  while  a  twist  is  applied  by  the  loading  mechanism  at  the 
other  end  of  the  input  bar.  The  release  of  the  torque,  that  is  stored  between  the  loading 
end  of  the  input  bar  and  the  clamp,  by  sudden  release  of  the  clamp  will  propagate  a  shear 
stress  pulse  along  the  bar  toward  the  specimen,  thus  loading  it  in  shear  at  high  rate. 


•  WEERASOORIYA,  T  Unpublished  Research,  1989. 

10.  SHRIVASTAVA.  S.  C.,  JONAS,  J.  J„  and  CANOVA,  G.  J.  Mcch.  Phys.  Solids,  v.  30,  1982.  p.  75. 

11  HARDING.  J.  Double-Notch  Shear  Testing  and  Punch  Loading,  High  Strain  Rate  Testing.  Mechanical  Testing.  Metal  Handbook,  American 
Society  for  Metals,  v.  8,  p.  228. 


2 


Figure  1  shows  a  photograph  of  the  THBar  apparatus,  and  a  schematic  diagram  of  the 
THBar  is  shown  in  Figure  2.  The  THBar  consists  of  two  25.4  mm  diameter  7075-T6  alumi¬ 
num  bars,  input  and  output  bars,  2438  mm  in  length.  These  two  bars  are  mounted  on  top  of 
a  square  box  beam,  using  oilite  bronze  bearings.  The  bearing  stands  are  designed  in  such  a 
way  that  the  bearings  have  the  freedom  to  move  in  any  direction;  i.e.,  translate  along  and 
rotate  about  three  mutually  perpendicular  axes  of  which  one  axis  is  lying  along  the  axis  of 
the  bearing.  Once  the  bearings  are  aligned,  they  are  locked  in  that  position.  Alignment  of 
the  bearings  is  performed  with  the  help  of  a  helium-neon  laser  and  a  steel  bar,  25.4  mm  in 
diameter  and  75  mm  in  length,  which  has  the  two  end  faces  precisely  machined  perpendicular 
to  the  axis  of  the  bar;  a  mirror  is  mounted  on  one  end,  and  a  cross  (two  lines  that  are  per¬ 
pendicular  to  each  other  and  passing  through  the  center),  are  scribed  at  the  other  end  of  the 
alignment  bar. 


Figure  1.  The  Torsional  Split-Hopkinaon  Bar  Apparatus. 


At  one  end  of  each  aluminum  bar,  a  hexagonal  slot  is  machined  to  enable  the  attachment 
of  a  hexagonal-end  specimen.  Specimens  are  attached  to  the  bar  with  12  slotted  screws,  two 
on  each  face  of  the  hexagon  (Figure  3).  Specimen  geometry  is  discussed  in  a  later  section. 
When  flat-end  specimens  are  tested,  these  bars  are  replaced  with  two  flat-end  bars.  In  this 
case,  using  an  epoxy,  the  specimen  is  glued  to  these  flat  ends  of  the  bars. 


Figure  2  Schematic  of  the  complete  Torsional  Hopkinson  Bar  Assembly, 


Figure  3.  Specimen-ends  of  the  bers.  in  this  picture,  hexagonal  end 
specimen  is  attached  to  the  bers  via  the  hexagonal  slots  in  the  bars. 


Shear  Pulse  Generating  System 

Torque  Generating  Mechanism 

A  torque  is  applied  to  the  opposite  end  from  the  slotted-end  of  the  input  bar.  via  a 
pulley  that  is  attached  to  the  bar  by  a  friction  fit.  Using  a  hydraulic  hand  pump  and  a 
cylinder,  the  pulley  is  rotated  to  apply  a  twist  to  the  input  bar.  The  schematic  and  photo¬ 
graph  of  this  loading  arrangement  is  shown  in  Figures  4  and  1.  tespectivcly.  The  pulley  is 
attached  to  a  cross  bar  by  cables.  The  cross  bar  is  pushed  down  by  the  hydraulic  cylinder, 
thus  rotating  the  pulley  that  is  attached  to  the  end  of  the  input  bar. 

Clamp 

A  photograph  of  the  clamping  arrangement  is  shown  in  Figure  5,  and  a  schematic  dia¬ 
gram  of  the  clamp  is  shown  in  Figure  6.  The  cLmp  consists  of  two  plates  that  are  posi¬ 
tioned  on  either  side  of  the  bar.  The  tops  of  the  plates  are  connected  by  a  2024-T4 
aluminum  pin,  19.1  mm  in  diameter,  with  a  12.2  mm  diameter  notch  at  the  center.  Both 
clamping  plates  are  mounted  in  a  carriage  that  is  free  to  move  on  ball  bearings  perpendicular 
to  the  bar,  with  pins,  allowing  the  plates  to  freely  rotate  about  these  pins.  While  clamping, 
application  of  lateral  forces  that  would  create  bending  and  axial  pulses  are  avoided  by  allow¬ 
ing  the  free  movement  of  the  clamp  relative  to  the  bar.  At  both  ends  of  the  clamp,  to 
reduce  any  bending  waves  emanating  from  the  release  of  the  clamp,  the  input  bar  is  mounted 
on  two  additional  38  mm  long  oilite  bronze  bearings,  which  are  located  on  either  side  of  ihe 
clamp.  In  addition  to  allowing  the  clamping  plates  to  rotate  about  the  attaching  pins,  the  pin 
of  one  plate  is  free  to  move  toward  the  bar.  Clamping  is  achieved  by  pushing  this  plate 
toward  the  bar  at  this  pin  by  a  hydraulic  cylinder  attached  to  a  hand  pump.  Further  pump¬ 
ing  of  the  hydraulic  oil  to  the  cylinder  will  fracture  the  notched  2024-T4  aluminum  pin  at  the 
top  of  the  bar,  suddenly  releasing  the  clamp  and,  therefore,  the  stored  torque  between  the 
clamp  and  the  loading  end. 
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Figure  4.  Schematic  of  th*  loading  anangemant 


Figure  5.  The  clamping  aeaembly  and  the  loading  pulley 


Lagrangian  Diagram  of  the  System 

Langrangian  ot'  x-t  diagram  of  the  system  is  shown  in  Figure  7.  The  position  of  the 
shear  wave  front  as  a  function  of  time  can  be  inferred  from  this  wave  diagram.  As  seen 
from  the  diagram,  at  the  release  of  the  clamp,  half  of  the  torque  stored  between  the  clamp 
and  the  loading  pulley  will  travel  as  a  shear  wave,  at  the  shear  wave  speed  (Cs)  of  aluminum, 
toward  the  specimen.  Shear  pulse  of  the  same  magnitude  will  travel  at  the  same  speed  in 
the  other  direction  unloading  the  torque  that  is  stored  at  that  end.  This  unloading  pulse  will 
reflect  completely  at  the  loading  pulley  (large  pulley  will  act  as  a  rigid  boundary)  and  travel 
toward  the  specimen  releasing  the  torque  to  zero;  therefore,  the  total  length  of  the  shear 
pulse  that  will  twist  the  specimen  is  equal  to  the  time  taken  by  the  unloading  pulse  to  travel 
toward  the  loading  pulle>  and  back  to  the  clamp.  By  allowing  the  clamp  to  be  relocated,  the 
length  of  the  pulse  can  be  changed.  To  allow  the  clamp  assembly  to  be  moved  and  boiled 
easily,  without  drastically  changing  the  alignment,  the  whole  clamp  assembly  is  mounted  on  a 
long  table  with  guidinii  rails. 
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Stress,  Strain,  and  Strain  Rate  Computation 

Magnitude  of  stress,  strain,  and  strain  rate  of  the  specimen  can  be  computed  form  the 
measured  shear  strain  records  of  the  input  and  output  bars  (Kolsky  analysis  method1).  Shear 
strain  of  the  specimen  (ys)  is  equal  to  the  relative  rotation  of  the  ends  of  the  gage  section 
of  the  specimen  divided  by  its  length: 


Ys  = 


Ds  0i  -  Ds  <t>2 
2Ls 


(1) 


where  <px  and  <p2  are  angles  of  twist  at  the  interface  between  specimen  (Figure  8)  and  input 
and  output  bars,  respectively,  Ls  is  the  gage  length  of  the  specimen,  and  Ds  is  the  mean 
diameter  of  the  thin  wall  of  the  specimen.  The  magnitude  of  <f>2  can  be  determined  from  yT: 
the  transmitted  shear  strain  at  the  surface  of  the  output  elastic  bar: 


_D  3<t>l  D  3<{>2 

^  "  2  dx  ~  2Cs  ~dT 


(2) 


where  Cs 
bar,  and  G 
rangement  yields: 


=Vg/p  is  the  shear  velocity  in  the  output 
and  p  are  shear  modulus  and  the  density  of 


bar,  D  is  the  diameter  of  the  output 
the  output  bar,  respectively.  Rear- 


<te  = 


(3) 
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Figure  9.  Schematic  of  the  specimen  bar  interfaces. 


The  magnitude  of  <px  can  be  determined  from  the  values  of  shear  strains,  the  incident  (•/[) 
and  reflected  (yR)  waves  at  the  specimen/input  bar  interface,  on  the  surface  elastic  input  bar: 


D  3<t>i  D  5<b i 
YI-7R  -  2  ~  2Q  7T 


(•*) 


and  integrating 
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I 

Ol  ==&■  I  [7l(0  -  7R(0]dt  . 
u  tf 

By  substituting  the  derivatives  of  <p\  and  <pi  from  Equations  3  and  5  in  Equation  1, 


(5) 


(6) 


By  integrating  the  above  relationship  for  specimen  strain  rate,  the  strain  of  the  specimen  can 
be  determined  as  a  function  of  time.  At  a  uniform  state  of  strain  in  the  specimen,  the  trans¬ 
mitted  torque  at  the  output  bar/specimen  interface  (T2)  is  approximately  equal  to  the  torque 
at  the  input  bar/specimen  interface  (Ti);  i.e.,  the  transmitted  shear  strain  pulse  is  equal  to 
the  sum  of  the  incident  and  reflected  shear  strain  pulses  at  the  specimen.  Therefore,  yT  ~ 

7t  +  7r  and  Equation  6  reduces  to: 

2C  D  (7) 

VsO)  7R«) 


Integration  of  Equation  7,  strain  of  the  specimen  can  be  obtained: 


Ys(0  = 


dt 


She^r  stress  (ts)  is  given  by: 


*s(t)  = 


2TS 

itDs2  ts 


(8) 


(9) 


where  t,  is  wall  thickness  and  Ts  is  average  torque  in  the  specimen:  Ts  is  the  average  of  Tj 
and  T2,  torques  at  the  input  bar  interface  and  output  bar  interface  of  the  specimen,  respec¬ 
tively.  These  torques,  Ti  and  T2,  can  be  computed  from  the  shear  strains  in  the  elastic  input 
and  output  bars  as  follows: 


GrcD3 

Ti - jy-Cn+Kk) 


(10) 


and 


T2= 


CmD3 
~ 16  ^ 


(ID 


As  discussed  earlier,  when  the  specimen  is  deforming  uniformly;  i.e.,  yx  —  7i  +  7r>  from 
Equation  9  using  Equations  10  and  11: 
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Measurement  of  Shear  Pulses  in  the  Bars 


Shear  stress  pulses  in  the  input  and  output  bars  are  measured  using  a  strain  gage  bridge. 
The  strain  gage  bridges  are  attached  to  each  bar  at  the  mid-points  (at  the  locations  A  and  B 
in  Figure  7).  By  locating  A  and  B  gages  at  the  mid-points,  both  transmitted  and  reflected 
waves  will  reach  the  gage  locations  at  the  same  instant;  also,  reflected  and  incident  pulse  will 
be  separated  in  time  in  the  trace  captured  at  the  A  gage.  In  addition,  another  strain  gage 
bridge  (C  in  Figure  7)  is  attached  between  the  clamp  and  the  loading  pulley;  this  bridge  is 
provided  to  measure  the  magnitude  of  the  applied  (stored)  torque  to  the  bar.  All  the  bridges 
consist  of  four  strain  gages.  Each  bridge  consisting  of  two  pairs,  one  pair  attached  diame¬ 
trically  opposite  the  other  on  the  surface  of  the  bar.  The  gages  are  oriented  at  a  45°  angle 
to  the  axis  of  the  bars;  gages  arranged  this  way  will  measure  only  shear  strains  by  excluding 
both  axial  and  bending  strains. 

Each  bridge  is  powered  by  a  separate  DC  power  supply,  an  HP6218C  constant  voltage 
(50V),  or  current  (0.2A)  power  supply,  in  which  the  voltage  can  be  set  with  a  potentiometer. 
Output  of  the  bridges  are  directly  fed  into  an  800  kHz  Nicolet  (Model  4562)  digital  oscillo¬ 
scope.  Power  supply  voltage  output  to  strain  gage  bridges  are  adjusted  to  achieve  a  chosen 
calibration  factor  and  is  discussed  in  the  Calibration  and  Data  Reduction  Method  Section. 

To  measure  axial  and  bending  pulses  that  will  be  generated  by  the  clamp,  four  additional 
strain  gages  are  attached  to  the  bar,  90°  apart,  at  each  A  and  B  location.  Two  diametrically 
opposite  gages  are  connected  in  series  to  act  as  the  active  part  of  a  single  gage  bridge-amplifier 
setup  (2  MHz  bandwidth).  Axial  stresses  are  measured  by  this  bridge  configuration.  The 
other  two  gages  that  are  180°  apart  are  a  part  of  a  two  active  gage  bridge-amplifier  unit 
(2  MHz)  that  is  set  up  to  measure  the  bending  stress  pulses  of  the  bar;  in  this  case,  the  two 
active  gages  are  located  at  opposite  locations  of  the  bridge. 

Calibration  and  Data  Reduction  Method 

Figure  9  shows  typical  waveforms  captured  by  the  digital  oscilloscope  at  A  and  B  gages. 
The  trace  from  A  gage  contains  both  the  incident  (vi(t))  and  reflected  (vR(t))  traces  sepa¬ 
rated  from  each  other  in  time  space.  The  trace  from  B  gage  represents  the  transmitted  pulse 
(vT(t))  through  the  specimen.  In  these  traces,  vertical  axes  represent  a  voltage  output  from 
the  strain  gage  bridge  that  is  proportional  to  the  shear  strain  at  the  surface  of  the  bars.  For 
a  four  gage  bridge,  the  shear  strain  (y)  at  the  surface  of  the  bar  is  given  by  the  expression: 

V  =  pg  »  (13) 

where  v  is  the  output  voltage  of  the  bridge,  F  gage  factor  of  the  strain  gages,  and  E  is  the 
input  DC  voltage  of  the  bridge.  By  using  Equations  13  and  7,  shear  strain  rate  of  the 
specimen: 


Ys(t)  = 


2CS  D$  2vr(i) 
Ls  D  FEa 


(M) 


where  subscript  R  represents  the  reflected  pulse  and  A  represents  the  gages  at  A.  Similarly, 
Equation  12  representing  shear  stress  will  become: 


ts(t)  = 


GD3  2vr(t) 
8Ds2ts  FEb 


(15) 


where  subscript  T  represents  the  transmitted  pulse  through  the  specimen  and  B  represents  the 
gages  at  location  B.  By  selecting  a  value  for  Ea/b,  (A  or  B)  calibration  factors  (i.e„  for  a 
unit  output  bridge  voltage)  for  strain  rate  and  stress  can  be  chosen.  For  chosen  calibration 
factors,  Ea  and  Eb  can  be  calculated  from  Equations  14  and  15.  With  these  input  bridge 
voltages,  when  a  known  precision  calibration  resistor  (Rc)  is  shunted  across  one  gage  of  the 
bridge,  the  output  bridge  voltages  (vr^)  can  be  computed  by  the  expression: 


VR/T  = 


REa/b 

4Fa/b(R  +  Rc) 


(16) 


where  R  is  the  resistance  of  a  single  gage. 


Shaw  Sira—  Pula—  at  A  and  B  La— Ilona 


Figure  9.  Shear  Stress  Pule—  at  A  and  B  Gage  Locations. 


Before  each  test,  by  tuning  the  pot  that  varies  the  magnitude  of  the  power  supply 
voltage,  the  calibration  factor  is  set  for  each  bridge  to  output  the  calculated  voltage  for  the 
calibration  resistor. 
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ACQUISITION  AND  ANALYSIS  OF  STRESS  PULSES 

Traces  of  the  stress  pulses  in  both  input  and  output  bars  are  captured  at  the  A  gage 
(incident  and  reflected)  and  B  gage  (transmitted)  using  the  digital  oscilloscope.  The 
oscilloscope  is  triggered  by  the  trace  at  gage  A  off  the  incident  pulse:  a  threshold  crossing 
with  a  prescribed  delay  of  -150  jus. 

The  waveforms,  captured  by  the  digital  oscilloscope,  are  transferred  to  a  microcomputer 
(AT  compatible)  via  an  RS232  serial  port  using  a  computer  program  written  in  WaveForm 
Basic  language  (see  Appendix  B  for  the  listing  of  the  program). 

Shifting  of  the  waveforms  in  time  space,  computation  of  strain  rate  (Equation  7),  stress 
(Equation  12),  strain  (by  numerically  integrating  strain  rate)  waveforms  from  captured  data  at 
A  and  B,  and  the  plotting  of  the  stress-strain  curves  are  performed  using  a  computer  program 
(written  in  TBASIC  language).  The  listing  of  the  computer  program  is  given  in  Appendix  B. 

EXPERIMENTS 

In  this  section  of  the  report,  preliminary  experiments  that  have  been  conducted  with  the 
THBar  are  discussed.  Also,  axial  and  bending  pulses  that  are  generated  by  the  clamp  are  pre¬ 
sented  and  discussed  in  this  section  of  the  report. 

Material 

Test  specimens  were  machined  from  oxygen-free-high-conductivity  (OFHC)  copper  (99.99% 
copper).  OFHC  material  was  received  as  one  inch  diameter  rods  in  the  work  hardened  condi¬ 
tion.  Before  machining  the  specimens,  copper  rods  were  annealed  at  400°C  vacuum  for  one 
hour.  This  annealing  heat  treatment  produced  grains  of  an  average  size  of  45  fi. 

Specimen  Geometry 

Geometry  of  the  test  specimen  is  shown  in  Figure  10.  Gage  section  of  the  test  specimen 
is  a  thin  wall  tube  (0.38  mm  wall  thickness)  of  0.254  mm  gage  length  and  of  outside  (Dso) 
and  inside  (Dsi)  diameters  of  10.16  mm  and  9.40  mm,  respectively.  A  thin  wall  tube  is  sand¬ 
wiched  between  two  hexagonal  flanges  which  are  used  to  attach  the  specimen  to  the  elastic 
input  and  output  bars  of  the  test  system.  Before  testing  the  specimens,  a  fine  line,  parallel 
to  the  axis  of  the  specimen,  is  scribed  on  the  inside  wall  of  the  specimen. 

The  geometry  of  the  specimen  is  identical  to  the  one  used  at  Brown  University  by 
Duffy.9  In  this  short  gage  length  specimen,  an  almost  homogeneous  state  of  strain  is 
achieved  after  a  few  reflections  of  the  loading  shear  stress  pulse.  In  this  specimen  geometry, 
although  the  plastic  zone  starts  at  the  flange-gage  section  interface,  the  plastic  zone  is  con¬ 
tained  until  it  spreads  gradually  through  the  specimen  and  engulfs  the  whole  gage  section.12 


12.  LEUNG,  E.  K.  C.  An  Elastic-Plastic  Stress  Analysis  of  the  Specimen  Used  in  the  Torsional  Kolsky  Bar.  J.  Appl.  Mech.,  v.  47,  1980,  p.  278. 
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DIM  INCHES 

A  0.600 
B  0.250 
C  0.100 

0  0.400±0.0005  DIA 
E  0.370±0.0005  DIA 


UNSPECIFIED  TOLERANCES  -  ±0.005 


Figure  10.  T «at  specimen  for  Torsional  Hopkinson  Bar. 


Tests 

Three  OFHC  specimens  were  tested  at  three  different  strain  rates,  400,  800,  and  1200  s'1, 
using  the  THBar.  The  specimens  were  attached  to  the  slotted  end  of  the  bars  by  screws. 

The  tests  were  initiated,  as  discussed  earlier,  by  first  clamping  the  bar  and  then  twisting  it  by 
the  loading  assembly  using  the  hydraulic  hand  pump  at  the  loading  arrangement.  While  clam¬ 
ping  the  bar,  it  was  made  certain  that  the  clamp  would  apply  its  two  clamping  forces  symme¬ 
trically  opposing  each  other,  thus  cancelling  any  lateral  force  components  that  would  apply 
any  bending  to  the  bar.  The  amount  of  twisting  was  indicated  by  a  digital  mili-voltmeter  that 
is  attached  to  the  strain  gage  bridge  at  the  location  C  (C  gage).  The  clamp  was  released  by 
further  hand  pumping  of  the  hydraulic  cylinder  at  the  clamp.  Shear  stress  pulses  at  the  A 
and  B  gages  were  captured  by  the  digital  oscilloscope. 

RESULTS  AND  DISCUSSION 

Captured  shear  pulse  traces,  at  A  and  B  strain  gage  locations,  for  one  of  the  tests  is 
shown  in  Figure  9.  As  shown  in  the  figure,  this  instrument  produces  an  approximately  con¬ 
stant  incident  pulse.  Rise  time  of  the  incident  shear  wave  generated  by  the  clamp  is  about 
45  fis. 
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Axial  and  bending  strain  pulses  at  the  A  gage  are  given  with  the  shear  strain  pulse  in 
Figure  11a  for  the  same  test.  Axial  and  bending  strain  pulses  that  are  generated  by  the 
clamp  are  appreciably  smaller  than  the  shear  strain  pulse  (less  than  10%). 

For  the  same  test,  after  transforming  measured  shear,  axial,  and  bending  strain  pulses  at 
the  B  gage  to  corresponding  stresses  through  the  specimen,  are  shown  in  Figure  lib.  In  the 
above  conversion  of  axial  and  bending  strains  at  B  gage  to  stresses,  the  following  expressions 
are  used: 


ED2 


and 


ED3  D,0 

Cb(t)  =  (bso^l5^£bT<t) 


(17) 


(18) 


where  aa  and  ab  are  axial  and  bending  stresses  at  the  specimen,  respectively,  eax  and  ebx  are 
axial  and  bending  strains  transmitted  to  the  output  bar,  respectively,  and  E  is  the  Young’s 
modulus  of  the  bar.  Bending  and  axial  stress  peaks  that  are  transmitted  through  the  specimen 
are  less  than  the  static  yield  of  the  annealed  OFHC  copper.  Moreover,  the  bending  stress 
pulse  travels  through  the  specimen  after  the  end  of  the  shear  pulse;  flexural  wave  speed  for 
the  aluminum  bar  supported  at  locations  L  (**  1300  mm)  apart  is  approximately  equal  to  2.9 
D/L  0.06)  times  the  shear  wave  speed  of  the  bar  for  the  zeroth  overtone  (free  end  condi¬ 
tion  is  assumed).  However,  the  axial  pulse  travels  faster  (*=  1.6  times)  than  the  shear  pulse 
and,  therefore,  it  reaches  the  specimen  earlier. 

The  captured  stess  pulses  at  gages  A  and  B  for  the  tests  are  transferred  to  the  microcom¬ 
puter.  The  reduction  of  these  data  are  accomplished  using  the  computer  program  in  Appendix  B. 
The  resulting  output  from  the  program,  stress-strain  plots  at  the  three  different  strain  rates, 
400,  800,  and  1200  s'1  are  shown  in  Figure  12.  Computed  shear  strain  rate,  stress  and  strain 
as  a  function  of  time  for  the  three  tests  are  shown  in  Figures  13,  14,  and  15,  respectively. 

Shear  stress-strain  curve  of  the  OFHC  copper  at  the  quasi-static  shear  strain  rate  of 
1  x  10‘3  s'1  is  also  shown  in  Figure  13.*  Stess-strain  curves  for  copper  at  strain  rates  of  400, 
800,  and  1200  s1  shows  a  strain  rate  sensitivity  of  this  material  in  the  range  of  rates  that  arc 
used  in  the  tests.  At  higher  strain  levels  («  35%),  under  high  strain  rate  conditions,  stress 
levels  seem  to  be  independent  of  the  strain  rates.  This  may  be  due  to  the  thermal  softening 
of  the  material  at  higher  strains  and  strain  rates.  However,  no  deformation  shear  bands  were 
present  in  the  gage  section.  High  strain  rate,  shear  stress-strain  data  (below  35%  strain  level) 
correlate  favorably  with  the  trend  from  the  quasi-static  rate  data  from  a  servohydraulic  test 
machine. 

It  is  noted  from  Figure  13,  after  a  60  pis  time  interval,  strain  rate  data  as  a  function  of 
time  for  the  three  tests  show  that  the  strain  rate  remained  reasonably  constant  during  the 
test.  This  60  pis  time  interval  corresponds  to  a  strain  of  approximately  1  to  5%,  depending 
upon  the  strain  rate  (see  Figure  15). 


•WEERASOORIYA.  T.  Unpublished  Research,  1989. 
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Figure  1 1  a.  Shear,  Axial,  and  Bending  Pulse*  at  the  A  Gage  Location. 
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Figure  1 1  b.  Shear,  Axial,  and  Bending  stresses  through  the  specimen 
from  captured  strain  pulses  at  the  B  Gage  Location. 
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Figure  1 3.  Strain  rate  as  a  function  time  for  the  three  tests. 
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Figure  1 4.  Stresses  as  a  function  of  time  for  the  three  tests. 
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Figure  1 5.  Strain  as  a  function  of  time  for  the  three  tests. 
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The  maximum  shear  strain  rate  that  can  be  achieved  in  this  THBar  configuration,  with 
the  specimen  of  gage  length  of  2.54  mm  and  the  present  clamp,  is  1400  s'1.  This  maximum 
limit  is  due  to  the  inability  to  apply  a  higher  clamping  force  to  prevent  the  slipping  of  the 
bar  at  the  clamp  for  higher  input  torques;  i.e.,  fracture  load  of  the  present  aluminum  fracture 
pin  is  low.  In  future  tests,  strain  rates  up  to  2500  s1  can  be  reached  with  the  present  specimen 
geometry  with  a  new  clamp  that  is  being  designed. 

CONCLUSIONS 

A  THBar  apparatus  was  developed  and  constructed  at  MTL.  This  apparatus  can  be  used 
to  generate  constant  shear  stress-strain  data  at  high  strain  rates,  up  to  1300  s'1,  for  the  deve¬ 
lopment  and  evaluation  of  constitutive  models  that  will  improve  the  description  of  the  shear 
behavior  of  the  materials.  In  addition,  the  shear  stress-strain  data  generated  by  this  apparatus 
can  be  compared  with  the  data  generated  by  the  compression  or  tension  Hopkinson  bars, 
using  the  equivalent  stress-strain  relationships  such  as  von  Mises.  These  comparisons  will 
allow  the  evaluation  of  existing  equivalent  stress-strain  relationships,  as  well  as  the  development 
of  better  ones. 

The  effect  from  the  bending  and  axial  pulses  generated  by  the  release  of  the  clamp  is 
determined  to  be  negligible.  Strain  rates  remain  reasonably  constant  after  60  n s  (i.e.,  5% 
strain  at  the  highest  achievable  strain  rate). 

The  apparatus  is  being  improved  to  have  additional  capabilities.  A  motor  drive  is  being 
developed  and  constructed  to  rotate  the  output  bar  at  the  far  end,  therefore,  allowing  quasi-static 
tests,  as  well  as  strain  rate  jump  (sudden  jump  of  strain  rate  from  the  quasi-static  to  a  higher 
rate  in  the  same  or  opposite  direction)  tests  to  be  conducted.  A  new  clamp  is  being 
developed  to  achieve  strain  rates  up  to  2500  s'1.  In  addition,  shear  bands  can  be  generated 
in  some  materials  with  this  new  clamp.13 
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APPENDIX  A.  COMPUTER  PROGRAM  TO  TRANSFER  WAVEFORMS 
FROM  DIGITAL  OSCILLOSCOPE  TO  COMPUTER 


10  TRANSFER  (D@,"a2”) 

20  REM  SELECT  FIRST  AND  LAST  POINTS  BY  OPERATOR-FOR  A-GAGE  DATA 
30  SETFXP 
40  SFTT  YP 

50  AGAGE@  =  CUT(D@) 

60  TRANSFER  (D@,”b2”) 

65  REM  USE  AGAGE  WINDOW  TO  CUT  B-GAGE  DATA 
70  BGAGE@  =  CUT  (D@) 

75  REM  PLOT  SELECTED  A-GAGE  AND  B-GAGE  DATA 
80  PLOT  AGAGE<2> 

90  PLOT  BGAGE@ 

95  REM  STORE  DATA  IN  AN  ASCII  FILE 
100  ASTORE  AGAGE@,BGAGE@,’TBARDATA” 

110  CLEAR  ALL 
120  CLOSE 
130  END 
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APPENDIX  B.  DATA  REDUCTION  PROGRAM 


1  !  HOPKINSON  TORSIONAL  BAR  DATA  ANALYSIS  PROGRAM 

2  !  WRITTEN  BY  TUSIT  WEERASOORIYA  (November  1989) 

3  SET  KEY  ON  I  GOTO  100 

4  PRINT  "Plot  STRESS  STRAIN  CURVE  on  ”,GinS 

5  SET  VIEWPORT  0,130,0,100 1  INPUT  PROMPT  "NEW  PAPER  AND  RETURN":AS 

6  Xmax  =  0.6  I  Xmin  =  0 1  Ymix  =  250 1  Ymin  =  0 

7  CLEAR  I  CAILStresa  strain  ILF  GinS="L"  THEN  COPY  I  END 

8  PRINT  "Plot  EPDOT  TIME  on  ".GinS 

1 0  SET  VIEWPORT  0,130,0.100 1  INPUT  PROMPT  "NEW  PAPER  AND  RETURN":A$ 

1 1  CLEAR  I  CALL  Epdotj  I  IF  GinS="L"  THEN  COPY  I  END 

1 2  PRINT  "Plot  STRAIN  TIME  on  ".GinS 

1 3  SET  VIEWPORT  0,130,0.100 1  INPUT  PROMPT  "NEW  PAPER  AND  RETURN“:A$ 

14  CLEAR  I  CALL  Epsilon_t  I  CALL  X_ticlab(l.E+06,l,"Time(uSec)") 

1 5  CALL  Y_ticlab(  100,1 , "Sir  ain(%)")  I  IF  GinS="L"  THEN  COPY  I  END 

1 6  PRINT  "Plot  STRESS_TIME  on_",Gin$ 

1 7  SET  VIEWPORT  0,130,0,100 1  INPUT  PROMPT  "NEW  PAPER  AND  RETURN":AS 
1  8  CLEAR  I  CALL  Sigma_t  1  CALL  X_ticlab(  1  .E+06, 1  ,’Time(uSec)") 

1 9  CALL  Y_ticlab(l,  1,  "Seres  s(MPa)")  I  IF  GinS="L"  THEN  COPY  I  END 

20  AS  =  HS  &  ".MOD"  I  PRINT  "Saving  Stress  Strain  curve  data  in  ".AS,"..." 

21  GOTO4180I  END 

24  [NP  PRO  "Send  graphics  to  the  Plotter,  LaserWriter,  Screen  (P/L/S) :  ":GinS 

25  GdeviceS  =  I  IF  GinS="P"  THEN  GdeviceS  =  “HPPLT=COMl:" 

26  IF  GinS="L"  THEN  GdeviceS  =  TBMPRN=D:M  I  SET  GRAPH  DEVICE  GdeviceS 

27  SET  TEXT  SIZE  2/2 1  SET  LINE  COLOR  1 1  CALL  Menu  I  END 

28  IF  GinS=  'L"  THEN  COPY  I  END 

40  CLEAR  I  CALL  Menu 

41  END 

100  ! - 

105  INITI  CLEAR 

110  US  =  CHRS(10)  I  AshiftS  =  "Y"  I  GinS  =  "S“  I  GdeviceS  *  "" 

120  PRINT  US 

210  OPEN  #6:"LPT1","W“ 

220  GinS  =  "S“  I INP  PRO  "Send  graphics  to  the  Plotter/Laser/Screen  (P/L/S)  :  ":GinS 
225  SELECT  CASE  GinS 
230  CASE'P" 

232  GdeviceS  =  "HPPLT=COMl:" 

233  CASE  “L" 

234  GdeviceS  =  "IBMPRN=D:" 

235  CASE  ELSE 

236  GdeviceS  =  "" 

237  END  SELECT 

238  SET  GRAPH  DEVICE  "" 

240  SET  TEXT  SIZE  2,2 1  SET  LINE  COLOR  1 

250  INPUT  PROMPT  "Enter  data  from  file  or  Default  file  (F/D) :  ":AS 

260  IF  AS="F  THEN 

265  PRINT  "Enter  the  Nicolel  Data  File  Name  (ex:  d:\fhame)  " 

270  INPUT  PROMPT "  (DO  NOT  ENTER  EXTENSION -assumed  PRN) :  ":HS 
280  ELSE 

290  HS  =  TB003" 

300  END  IF 

310  HeS  =  HS  A  "  PRN" 

330  Dbar  =  1.0  I  Dout  =  0.4  I  Din  =  0.37  !  Lspec  =  0.1  !  inches 
335  Cl  =  5036  I  Cs  =  3100  !  m/sec. 

340  OPEN  #7:HeS,"F" 

350  AS  =  "K"!INP  PRO  "Bar  and  Spec,  dimensions  from  file  or  key  board -inches  (F/K) :  ":AS 

360  IF  AS="F  THEN 

370  INPUT  #7:Dout,Di, Lspec 

380  ELSE  IF  AS="K" THEN 
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390  ! INPUT  PROMPT  "Enter  Dout,  Din  and  Spec,  length  (in.):  ":Dout,Din,Lspec 

400  ELSE 

410  GOTO  350 

420  END  IF 

425  Dm  spec  =  (Dout+Din)/2.0  I  Th  =  (Dout-Din)/2.0  !  Th  =  Th*0.0254 
430  Dbar  =  Dbar*0.0254  1  Dmspec  =  Dmspec*0.0254  I  Lspec  =  Lspec*0.0254 

435  Calib  a  =  25000*(0.1*Dmspec/(0.385*Lspec)) 

436  Calib.b  =  4.E+04*(0.385A2*0.015*0.0254A3/(DmspecA2*Th)) 

440  Rbar  =  Dbar/2  I  Areab  =  Pi/4*DbarA2  I  Areas  =  Pi/4*  Dmspec  A2 
460  Delb  =  48.5*0.0254  I  Dele  =  485*0.0254  I  Dels  =  59*0.0254 
465  E  =  7.1E+10  I  G  =  2.69E+10  !  MPa 

470  Ro  =  2800  I  Adelt  =  5.E-05 

475  Cl  =  SQR(E/Ro)  I  Cs  =  3100  1  CS  =SQR(G/Ro)  -  m/sec. 

477  Tshift  =  2*Delb/Cs 
480  [read  =  4094 

490  1  READ  NICOLET  TWO  WAVE-FORMS - 

495  DIM  Temp(2,Iread]  I  Temp  =  0 

507  ON  EOF(7)  GOTO  590 

510  INPUT  #7 :T0.T emp(  1 , 1  ] ,T0,Temp( 2, 1  ] 

520  Xmin  =  T0/1.E+O6  !  Seconds 

530  Ofsetl  =  Temp(l,ll  I  Ofset2  =  Temp(2,l) 

545  INPUT  #7:T,Temp[  1 ,2].T,T emp(2.2] 

546  Delt  =  (T-T0)*  1 JE-06  !  Seconds 
560  FOR  1  =  3  TO  Iread 

570  INPUT  #7:T,Temp{U].T.Temp[2,I] 

580  NEXT  I 
590  CLOSE 7 
595  Imax  =  1-1 

597  DIM  Wform(2Jmax],TempO[Imax]  1  Wform  =  0  I  Tempo  =  0 

600  FOR  1=1  TO  Imax 

610  Wfonn[lJ]  =  Tempi  1J]  -Ofsetl 

620  Wfoim[2J]  =  Temp(2.I]-Ofset2 

630  Temp0(I]  =  Wform(l,I] 

660  NEXT  I 

661  DELETE  Temp 

662  DIM  Ep{3 Jmax]  I  Ep  =  0 

670  ! - 

680  Xmax  =  T/l .E+06-Xmin  I  Xmin  =  0  ITxmax  =  Xmax  ITxmin  =  Xmin 

690  Xran  =  Xmax-Xmin 

700  CLEAR 

760  CALL  D_windows 

780  Delt  =  Xrart/(lmax-l )  I  Tmax  =  Txmax  I  Imaxp  =  INT(Tmax/Deli) 

790  CALL  P  MIN(TempO,Ymin,Ypos)  I  Ymax  =  -Ymin  !  CALL  P_MAX(TempO,Ymax,Ypos) 
800  IMAGE  2D.6(4D.4D) 

810  Delx  =  ABS(Xmax-Xmin)  I  Dely  =  ABS(Ymax-Ymin) 

820  SET  WINDOW  Xmin-0.2*Delx.Xmax+0.1*DeIx.Ymin-0.1*Dely.Ymax+0.1  *Dcly 

830  !  DRAW  TWO  NICOLET  WAVEFORMS  IN  ACTUAL  TIME  SPACE - 

840  SET  VIEWPORT  0.6550.100 1  SET  LINE  COLOR  1 1  SET  LINE  STYLE  0 
850  MOVEGDU  129.2 1  SET  TEXT  ALIGN  55  I  SET  TEXT  SIZE  15,15 
860  PLOT  TEXT  GDU  HS I  SET  TEXT  ALIGN  0,0 1  SET  TEXT  SIZE  12 
no  CALL  X  ticlab(  1  .E+06, 1  .'Time(uSec)")  1  CALL  Y_ticlab(100.1."Stram(%)") 

880  AXIS  Xinc,Yinc,0.0,3.3 
890  FOR  J  =  I  TO  2 

900  SET  LINE  COLOR  J+l  I  SET  TEXT  COLOR  J+l 

910  FOR  1  =  l  TO  Imax 
920  X  =  Xmin+(I-l)*Delt 

930  Y  =  WformfJ.I] 

940  IF  1=1  THEN  MOVE  X.Y  ELSE  DRAW  X.Y 

950  NEXT  I 

960  IF  J=1  THEN  TiileS  =  “Ea," 

970  IF  J=2 THEN  TitleS  =  "  Eb“ 

990  MOVE  GDU  30+J*  12,90 

1000  PLOT  TEXT  GDU  TitleS 

1010  NEXTJ 
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1015  !  SHIFT  REFLECTED  AND  TRANSMITTED  WAVEFORMS  AT  B  BY  TRAVEL  TIME 
1017  MOVE  Tshift.Ymin/2  I  DRAW  Tshift,Ymax/2 
1020  IF  AshiftS='Y"THEN 

1022  PRINT  'Twice  the  Time  (uS)  from  A  ->  Specimen  ";Tshift*l.E+06; 

1 024  INPUT  PROMPT  "Change  the  time  (Y/N)  ?":AshiftS 
1026  IF  Ashift$="Y"  THEN 

1 028  INPUT  PROMPT  "Enter  Time  (uS)  =  ":Tshift  I  Tshift  =  Tshift/1  .E+06 

1030  ELSE 

1032  SET  GRAPH  DEVICE  GdeviceS 

1034  END  IF 

1039  GOTO  700 

1040  END  IF 

1042  Ishift  =  INT(Tshift/Delt) 

1045  FOR  I  =  Ishift  TO  Imax 
1050  J  =  I-Ishift-t-1 

1060  Ep[2,J]  =  WformflJ] 

1065  Ep{3,Jl  *  Wform(2.I] 

1070  NEXTI 

1072  FOR  1=1  TO  Imax 

1073  EpIU]  =  Wform(lJ] 

1074  NEXTI 

1080  Itest  =  Imax-Ishift+l 

1110  !  PLOT  THE  SHIFTED  WAVEFORMS - 

1120  SET  VIEWPORT  0,65,0.50 1  SET  LINE  COLOR  1 1  SET  LINE  STYLE  0 
1130  MOVE  GDU  129,2 1  SET  TEXT  ALIGN  5,3  I  SET  TEXT  SIZE  1.5,1 .5 
1140  PLOT  TEXT  GDU  HS I  SET  TEXT  ALIGN  0,0 1  SET  TEXT  SIZE  2£ 

1 150  CALL  X_ticlab(l .E+06, l.'Time(uSec)")  I  CALL  Y  ticlab(100.1."Strain(%)") 

1160  AXIS  Xinc, Yinc,0,0.3 ,3 
1170  FOR  J  =  1  TO  3 

1180  SET  LINE  COLOR  J+l  I  SET  TEXT  COLOR  J+l 

1190  FOR  1=1  TO  Itest 
1200  X  =  Xmin+(I-l)*Delt 

1210  Y  =  EpfJJ] 

1220  IF  1=1  THEN  MOVE  X,Y  ELSE  DRAW  X,Y 

1230  NEXTI 

1 240  IF  J=  1  THEN  TitleS  =  "  E-ins," 

1 250  IF  J=2  THEN  TitleS  =  "  E-ref," 

1 260  IF  1=3  THEN  TitleS  =  "  E-tm" 

1 270  MOVE  GDU  30+J*1Z90 1  PLOTTEXT  GDU  TitleS 

1280  NEXTJ 

2330  '.Eplns,  EpRef,  EpTm  and  (EpIns+EpRef) - 

2340  SET  WINDOW  Xmin-0.2*Delx,Xmax+0.1*Delx,Ymin-0.1*Dely,Ymax+0.1*Dely 

2350  SET  VIEWPORT  65,130,0,50 1  SET  LINE  COLOR  1 1  SET  LINE  STYLE  0 

2360  MOVE  GDU  129,2 1  SET  TEXT  ALIGN  53  I  SET  TEXT  SEE  1.5,13 

2370  PLOT  TEXT  GDU  HS  I  SET  TEXT  ALIGN  0,0 1  SET  TEXT  SIZE  Z2 

2380  CALL  X  ticlab(l.E+06,l,'Time(uSec)")  (  CALL  Y  ticIab(I00,l.“Strain(%)") 

2390  PLOT  AXIS  Xinc,Yinc,0,0,3,3 
2400  FOR  J  =  1  TO  4 

2410  SET  LINE  STYLE  01  SET  LINE  COLOR  J+l  I SETTEXT  COLOR  J+l 
2420  FOR  I  =  I  TO  Itest 
2430  X=Xmin+a-l)*Delt 

2440  IF  Jot  THEN  Y  =  EptfJ]  ELSE  Y  =  -(Ep{l,n+Ep(ZIl) 

2450  IF Iol  THEN  DRAW  X,Y  ELSE  MOVE  X.Y 

2460  NEXTI 

2470  IF  J=  1  THEN  TitleS  =  "  Eplns," 

2480  IF  J=2  THEN  TitleS  =  "  EpRef," 

2490  IF  J=3  THEN  TitleS  =  "  EpTm,” 

2500  IF  J=4  THEN  TitleS  =  "  EpIns+EpReF' 

2510  MOVE  GDU  5+<J-l)*23.90 1  PLOT  TEXT  GDU  TideS 

2520  NEXTJ 

2525  IF  GinS="L"  THEN  COPY  else  INPUT  PROMPT  New  paper  AND  Press  ENTER  "  AS 
2540  CLEAR  I  CALL  D_ windows 

2550  !  Computation  of  Sigma,  Epdot  AND  Epsilon . 

2555  DELETE  Wform.TempO 
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2S60  DIM  Sigma[Itest],Epdot(Itest],Epsilon[Itest],Asigma[Itest] 

2570  FOR  I  =  1  TO  best 

2580  Asigmad]  =  -Calib  b*(-(Ep(l,I]+Ep{2,I])+Ep[3.I])/2 

2590  Sigmall]  =  -Ep[3.I]*Calib  b  !  ((G*dbar)  /  ( (8*dmspec*2)  *th) ) 

2600  Epdot(I]  =  -Calib  a*Ep{2J]  !  ((2*Cs*dmspec)/(lspec*dbar)) 

2610  NEXTI 

2620  !  Compute  Epsilon  by  Integrating  Epdot - 

2630  Epsilon{l]  =  (Epdot(l]+Epdot(2])*04*Delt 
2640  FOR  1  =  2  TO  Itest-1 

2650  Epsilonfl]  =  Epsilonfl- ll+Delt*04*(Epdot[T]+Epdot[I+l]) 

2660  NEXTI 

2670  Epsilonfltest]  =  Epsilon(Itest-l] 

2680  !  Epdot-Time  Plot - 

2690  SET  VIEWPORT  0,6540. 100 1  CALL  Epdot_t 

2700  !  Sigma-Time  Plot - 

2710  SET  VIEWPORT  0,65,50.100  I  CALL  Sigma.t 

2720  SET  VIEWPORT  65,13040,100  I  CALL  X_ticIab(l.E+06,l,Time(uSec)") 

2730  CALL  Y_ticlab(l.l,,'Stress(MPa)'')  I  CALL  Sigmaj 

2740  !  Epsilon-Time  Plot - 

2750  SET  VIEWPORT  0,65.50.100  I  CALL  Epsilon.t 

2760  SET  VIEWPORT  65.130,040  I  CALL  X_ticIab(l.E+06, l.Time(uSec)") 

2770  CALL  Y_ticlab(100,l,"Sirain(%)H)  I  CALL  Epsilon_t 

2780  !  Sigma-Epsilon  plot - 

2790  SET  VIEWPORT  0,65,0.50 

2800  CALL  P_MIN(Sigma.Ymin,Ypos)  I  CALL  P_MAX(Sigma,Ym«x,Ypos) 

2810  CALL  P_MIN(Epsilon,Xmin,Ypos)  I  CALL  P_MAX(Epsilon,Xmax,Ypos) 

2825  CALL  Stress_strain 

2826  IF  GinS=”LH  THEN  COPY  I  PRINT  LfS  I  PRINT  "F10  for  Menu" 

2827  END 

2828  !  Sigma-Epsilon  plot - 

2829  SUB  S  tress  ^strain 

2830  Debt  =  ABS<Xmax-Xmin)  I  Dely  =  ABS(Ymax-Ymin) 

2840  SET  WINDOW  Xmin-0.2»DclxJCm*x+0.1*Delx,Ymm-0.2*Dely,Ymax+0.1*Dely 

2850  SET  LINE  STYLE  01  SET  LINE  COLOR  1 

2860  CALL  X  Uclab(100.1,"Strain(%n  1  CALL  Y  ticlab(l.l.,  Stress(MPa)“) 

2870  AXI  Xinc.Yinc.0,044  I  SET  UN  STY  0 1  SET  UN  COL  5 1  SET  TEX  COL  5 

2880  FOR  1=1  TO  Itest-1 

2890  X  =  Epsilon! I]  I  Y  =  Sigmafl] 

2900  IF  Iol  THEN  DRAW  X.Y  ELSE  MOVE  X.Y 

2910  NEXTI 

2920  !  SET  LINE  STYLE  II  SET  UNE  COLOR  1 

2930  !  FOR  I  =  1  TO  itest-1 

2940  !  X  =  Epsilon[f]  I  Y  =  Asigma(T) 

2950  !  IF  Iol  THEN  DRAW  X.Y  ELSE  MOVE  X.Y 

2960  (NEXTI 

2970  TideS  =  ”Sigma(MPa)  Epsilon(%)"  I  SET  TEXT  COLOR  1 
2980  SET  UNE  STYLE  01  MOVE  GDU  40,90 1  PLOT  TEXT  GDUTtdeS 

2990  MOVE  GDU  129 2 1  SET  TEXT  ALIGN  54  I  SET  TEXT  SIZE  14,14 
3000  PLOT  TEXT  GDU  HSISET  TEXT  ALIGN  0,0 1  SET  TEXT  SIZE  24 

3020  END  SUB 

3690  I - 

3700  SUB  D  windows 

3710  CLEAR  I  SET  UNE  STYLE  01  SET  UNE  COLOR  2 
3720  SET  VIEWPORT  0,130,0,100  I  SET  WINDOW  0,130,0,100 

3730  MOVE  0,0  I  DRAW  0.0;130,0;130,100;0,100;0,0 

3740  MOVE  65,0  I  DRAW  65, 100  I  MOVE  0,50  I  DRAW  13040 

3750  END  SUB 

3760  ! - - 

3770  SUB  X  ticlablXmul.Tcol.XlabelS) 

3780  SET  TEXT  COLOR  Tcol  I  SET  CUP  OFF 

3790  X  =  LOG10(Xmax)  I  Xinc  =  1NT(X)  I  Xinc  =  lO^Xinc 

3800  IF  Xmax>5*Xinc  THEN 

3810  Xiinc  =  2 

3  820  ELSE  IF  Xmax>24*Xinc  THEN 
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3830  Xiinc  =  1 

3840  ELSE 

3850  Xiinc  =  0.5 

3860  END  IF 

3870  1=01  SET  TEXT  ALIGN  3,0 

3880  DO 

3  890  MOVE  I*Xinc,-0.Q5*Dely  I  PLOT  TEXT  GDU  I*Xinc*Xmul 

3900  IF  (I+Xiinc)*Xin»Xmax  THEN  EXIT 

3910  I  =  I+Xiinc 

3920  LOOP 

3930  SET  TEX  AU  3,3  I  MOV  Xmin+0.5*Delx.-0.1*Dely  I  PLO  TEX  GDU  XlabelS 
3940  SET  TEXT  AUGN  0,0 1  SET  CLIP  ON 

3950  END  SUB 

3960  ! - 

3970  SUB  Y  ticlab<Ymul.Tcol,Ylabel$) 

3980  SET  TEXT  COLOR  Tcol  I  SET  CLIP  OFF 

3990  Y  =  LOG10(Ymax)  I  Yinc  =  INT(Y)  I  Yinc  =  10*Yinc 

4000  IF  Ymax>5*Yinc  THEN 

4010  Yiinc  =  2 

4020  ELSE  IF  Ymax>2.5*Yinc  THEN 

4030  Yiinc  =  1 

4040  ELSE 

4050  Yiinc  =  0.5 

4060  END  IF 

4070  I  =  01  SET  TEXT  AUGN  5.3 

4080  DO 

4090  MOVE  OJ'Yinc  I  PLOT  TEXT  GDU  I*Yinc*Ymul 

4100  IF  a+Yiinc)*  Yino  Ymax  THEN  EXIT 

4110  I  =  1+ Yiinc 

4120  LOOP 

4130  SET  TEXT  AUGN  3.3 1  SET  TEXT  ANGLE  90 

4 140  MOVE  -0.17*Delx,Ymin+0-5*Dely  I  PLOT  TEXT  GDU  YlabelS 

4150  SET  TEXT  ANGLE0I  SET  TEXT  ALIGN  0,01  SET  CUP  ON 

4160  END  SUB 

4170  ! - 

4180  AS  =  HS  A  ".MOD"  I  OPEN  #7:AS,"F" 

4190  FOR  1=1  TO  Itest 

4200  !  IF  Epsilon(H>0.6  OR  Sigma{Il>100  THEN  GOTO  4230 

4210  PRINT  #7:(i-l)*delt*le6;EpsiIon[I];Sigma(I];epdot[i] 

4220  NEJCTI 
4230  CLOSE 7 
4240  END 

4250  ! - 

4260  SUB  Epdoc  t 

4270  CALL  P_MIN(Epdoi,Ymin,Ypos)  I  CALL  P_MAX(Epdot,Ymax,Ypos) 

4280  Xmax  =  Txmax  I  Xmin  =  Txmin 

4290  Debt  =  ABS(Xmax-Xmin)  I  Dely  =  ABS(Ymax-Ymin) 

4300  SET  WINDOW  Xnin-0.2*Delx,Xmax+0.1*Delx,Ymin-0.2*Dely,Ymax+0.1*Dely 

4310  SET  UN  COL  2 1  SET  UN  STY  01  CALX  tklab(I.E+06,l."Time<uSec)") 

4320  CAL  Y_ticlab(l,2,"Strim  Rate(l/Sec)")  I  PLO  AXI  Xinc,Yinc,0,0,3,3  I  SET  TEX  COL  2 

4330  SET  UNE  STYLE  01  SET  LINE  COLOR  2 1  SET  TEXT  COLOR  2 

4340  FOR  I  =  1  TO  Itest 

4350  X  =  Xmin+(I-l)*Delt  I  Y  =  Epdotfl] 

4360  IF  Iol  THEN  DRAW  X.Y  ELSE  MOVE  X.Y 

4370  NEXT! 

4380  TitleS  =  "StiainRate  t“ 

4390  MOVE  GDU  87,90 1  PLOT  TEXT  GDU  TitleS 

4400  MOVE  GDU  129,2 1  SET  TEXT  ALIGN  5.3  I  SET  TEXT  SIZE  1.5. 1.5 
4410  PLOT  TEXT  GDU  HS !  SETTEXT  AUGN  0,0 1  SET  TEXT  SIZE  22 

4420  END  SUB 

4430  !  MENU . . 

4435  SUB  Menu 
4440  PRINT  US 

4450  PRINT  "FI  Plot  STRESS  STRAIN  Curve” 
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4460 

4470 

4480 

4490 

4500 

4505 

4510 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 


PRINT  "F2  -  Plot  EPDOT_TTME  Curve" 

PRINT  "F3  -  Plot  STRAIN  TIME  Curve" 

PRINT  "F4  -  Plot  STRESS  TIME  Curve" 

PRINT  "F5  -  Write  STRESS_STRAIN  Date  to  a  Disk  File" 
PRINT  "F6  -  Select  Graphic  Device" 

PRINT  "F7  -  Print  Graph  to  the  LaserWriter" 

PRINT  "F10  -  MENU" 

END  SUB 


SUB  Sigma_t 

CALL  P_MIN(Asigma,Yil,Ypos) !  CALL  P_MAX(Asigma,Yal,Ypos) 

CALL  P_MIN(Sigma,Yi2,Ypos)  I  CALL  P_MAX(Sigma,Ya2.Ypos) 

Ymin  =  MIN(Yi  1 , Yi2)  I  Ymax  =  MAX(Yal,Ya2) 

Debt  =  ABS(Xmax-Xmin)  1  Dely  =  ABS(Ymax-Ymin) 

SET  WINDOW  Xmin-0.2*Delx,Xmax+0.1*DeIx,Ymin-0.2*Dely,Ymax+0.  l*Dely 

SET  UNE  COLOR  3 1  SET  LINE  STYLE  0 

Y  =  LOG10(Ymax)  I  Yinc  =  INT(Y)  I  Yinc  =  10*Yinc 

X  =  LOGlOCXmax)  I  Xinc  =  INT(X)  I  Xinc  =  10*Xinc 

PLOT  AXIS  Xinc,  Yinc, 0,0,3 ,3  I  SET  LINE  COLOR  3  I  SET  TEXT  COLOR  3 

FOR  I  =  l  TO  Itest 

X  =  Xmin+(I-l)*Delt  I  Y  =  Sigmafl] 

IF  Iol  THEN  DRAW  X.Y  ELSE  MOVE  X.Y 
NEXT  I 

TideS  =  “Stress_Time"  I  SET  UN  STY  0 1  MOV  GDU  87,85  I  PLO  TEX  GDU  TideS 
SET  TEXT  COLOR  1 1  SET  LINE  COLOR  1 
FOR  I  =  1  TO  Itest 

X  =  Xmin+(I-l)*Delt  I  Y  =  Asigmafl] 

IF  Iol  THEN  DRAW  XY  ELSE  MOVE  X.Y 
NEXT  I 

TideS  =  "AvStress_Time"  I  MOVE  GDU  87.80 1  PLOT  TEXT  GDU  TideS 
MOVE  GDU  129,2 1  SET  TEXT  AUGN  53  I  SET  TEXT  SIZE  13.13 
PLOT  TEXT  GDU  HS I  SET  TEXT  AUGN  0,0 1  SET  TEXT  SIZE  22 
END  SUB 


SUB  Epsikm_l 

CALL  P_MIN(Epsilon.Ymin,Ypos)  i  CALL  P_MAX(Epsilon.Ymax.Ypos) 

Debt  =  ABS(X/nax-Xmin)  I  Dely  =  ABS( Ymax -Ymin) 

SET  WINDOW  Xmin-0.2*Delx,Xmax+0.1*Delx,Ymin-0.2*Dely,Ymax+0.1*Dely 

SET  UNE  COLOR  4 1  SET  UNE  STYLE  0 

PLOT  AXIS  Xinc,  1,0.033  I  SET  TEXT  COLOR  4 

FOR  I  =  1  TO  Itest 

X  =  Xmin+<I-l)*Delt  I  Y  =  EpsilonP] 

IF  Iol  THEN  DRAW  X.Y  ELSE  MOVE  X.Y 
NEXT  I 

TideS  =  "Strain  Tune"  I  MOVE  GDU  87,75  I  PLOT  TEXT  GDU  TideS 
MOVE  GDU  1293 1  SET  TEXT  AUGN  53  I  SET  TEXT  SIZE  13,13 
PLOT  TEXT  GDU  HS  I  SET  TEXT  AUGN  0.0 1  SET  TEXT  SIZE  22 
END  SUB 
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